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TECHNICAL NOTE NO. 1852 


REINFORCED CIRCULAR CUTOUTS IN PLANE SHEETS 
By H. Rei saner and M. Morduchow 


SUMMARY 


The problem treated here is to design the reinforcement of a cutout 
in a plane sheet in such a way that it is as nearly as possible equivalent 
to the part of the structure which has been cut out. A perfect equivalence 
would mean that the stresses and displacements of the structure remain the 
same as those which would have appeared without the cutout. 

General formulas are developed for the circumferential distribution of 
the cross-sectional moment of inertia I p and of the area Ay of a 

circular reinforcement required for perfect equivalence. These formulas 
are then applied to some cases of external edge tractions: Hydrostatic 

stress, pure shear, uniaxial tension, and pure bending. It is found that 
in the first two cases, the required cross sections are physically possible 
(i.e., I r and Ay come out positive), although the required moment of 
inertia is in some cases found to be quite high in comparison with the 
required area. 

In the cases of uniaxial tension and pure bending, it is shown that 
constraint stresses, that is, additional stresses in the sheet due to the 
reinforced cutout, are practically unavoidable. Simple formulas are 
developed for calculating these "constraint" stresses for any given (constant) 
cross-sectional characteristics of the reinforcement ring. These formulas 
are derived on the basis of the assumption that the constraints diminish 
sufficiently rapidly with radial distance from the cutout so as to have 
little effect at the external edges of the sheet. 

To check the influence which actual boundary conditions might have on 
the practical validity of these formulas, a test was made on a plane sheet 
with a reinforced circular cutout subjected to a tensile load causing 
constant displacements at the loaded edges. It was found that the values 
of the strains calculated from the exact formulas developed in this report 
for an infinite sheet were fairly similar to the values of the measured 
strains in the specimen, except along the loaded edges, where the actual 
strain decreased more rapidly toward the center of the sheet than the 
calculated strains. This discrepancy must be due, at least in part, to 
the actual condition of constant displacements at the loaded edges of 
finite length instead of, as in the analytical formulas, constant stress 
at the remote loaded edges of infinite length. 


2 


NACA TN No. 1852 


INTRODUCTION 


This investigation is concerned with the problem of the reinforcement 
of cutouts in plane sheets. Such cutouts may serve either to make a 
structure lighter or to provide space for personnel or accessories. The 
problem of designing a reinforcement ring around the cutout so that it will 
be elastically equivalent to the portion of the sheet removed is treated 
here by an inverse method. The treatment is given first in general terms 
and then applied in detail for several types of edge loading. 

A test was made on a plane sheet with a reinforced cutout to check 
the effect of some of the simplifying assumptions on the formulas developed. 

This investigation was carried out at the Polytechnic Institute of 
Brooklyn under the sponsorship and with the financial assistance of the 
National Advisory Committee for Aeronautics. Grateful acknowledgement is 
hereby also made to Professor N. J. Hoff and to Dr. B. Boley for giving 
the authors the benefit of the experience of the Laboratory for Aircraft 
Structures of the Polytechnic Institute of Brooklyn. 


SYMBOLS 


Ar 

a 

a i 

b 

c mn 

d 


e 

E, 


E, 


h 


effective area of reinforcement ring with rivet holes 

cross-sectional area of cutout reinforcement ring 

radius of center line of ring 

radius of outer circumference of ring 

half width of sheet 

arbitrary constants in stress function 
radial width of ring 

distance from edge to neutral fiber of reinforcement ring 

modulus of elasticity of ring material 

modulus of elasticity of sheet material 

bending moment in cross section of ring 

height of web 
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I 

moment of inertia; with subscripts e, F, and 

the effective moment of inertia, moment of inertia of 
the flanges, and moment of inertia of the web, 
respectively 

J r 

cross-sectional moment of inertia of ring 

i 

f E r T r 

dimensionless parameter of moment of inertia - y-2k(l + V 

\^s a-H 

l 

half length of sheet 

^mn 

abbreviations for products of c^ and powers of a-^ 

M 

bending moment about z-axis 

N 

radial shear force on cross section of reinforcement ring 

n 

order of terms of stress function 

P 

load 

P n ( r ), Qn(r) 

coefficients in stress function (cf. equation (3 a )) 

r 

radial distance from center of cutout 

r o 

radial distance from center of cutout to neutral fiber 

t 

thickness of sheet 

T 

normal stress resultant 

u, V 

displacements in x- and y-directions, respectively 

TV, ^ 

same displacement system in radial and transverse 
directions, respectively 

y, z 

Cartesian coordinates, measured from center of cutout 

Y, Z 

tangential and radial loads on ring per unit of 
circumference 

a 

( E Aj* \ 

dimensionless parameter of area of ring ^ ^ ^?2(l + ^ ) 

?rq) 

shear strain 


6 


dimensionless parameter of width of ring (d/2a) 
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K r 


V 


P 


°x ’ °y> T 


*y 


a 


r > 



T 

rep 


T 

O 


9 

* 


strain in x- and y-directions, respectively 
radial and transverse strain, respectively 
change of curvature of ring 
Poisson's ratio 
radius of curvature 

normal and shear stresses in sheet (Cartesian coordinates) 

radial, transverse, and shear stresses in sheet (polar 
coordinates) 

original shear stress 

angular polar coordinate (fig. 1) 

stress (Airy) function 


Sub- subscripts: 

0 original 

1 constraint 

T total 


GENERAL EQUATIONS AND GENERAL SOLUTIONS 


In order to analyze the effect of a circular cutout in a plane sheet, 
the stresses, displacements, and strains in such a sheet may he represented 
in a polar -coordinate system (fig. l) • The original stresses and displace- 
ments, before being disturbed by the cutout, will then appear in the form 

\ 


9 

^o 


1 o 
x o 


( 1 ) 
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The cutout edges are assumed to carry reinforcement strips or angles. It 
is the purpose of the analysis of this paper to choose the elastic properties 


of such an edge reinforcement so that the increments u r , 


and 


r rcf>L 


to the original stresses (the original stresses being given by the sheet 
edge loading) become zero or as small as possible. 


Stresses 

The stress in a plane sheet can be expressed by a stress function 
i(r, cp) such that '' 


o r = r~ 



where the primes and dots are defined by 




J 


( 2 ) 


d_ 

dr 

a 


I 


The function t, moreover, must be such that it satisfies the equation 


AA| r = 0 


where 


A\|r 


d 2 \|r a 2 t 

dx 2 dy 2 

1 8/ djjr\ 1 d 2 f 

r dr V dr/ r^ dq? 


( 3 ) 


The complete solution of equation ( 3 ) for the case considered here, 
in which the edge stresses taken around the circumference of the circular 
cutout have zero resultant, may be written in the form: 
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t = y 1 P n cos ncp + y^ sin 


ncp 


(3a) 


n=0 


n=l 


where 


\ 


P Q = c 10 + c 20 r + c 3Q log e r + c^r log e r 


P 1 = c ll r + c 21 r3 + c 3 i r " 1 + c hl r lQ g e r 


P = c, r 11 + c 0 r” n + c_ r n+2 + c, r“^ n “ 2 ^ 
n>2 In 2n 3n J+n 




(3h) 


J 


The expressions for the functions Q are the same, except with 
possibly different numerical values of the constants c. 

From equations (2) and (3a) it follows that 


°r ' X (-* 2 *- 2 *n + r ' lp n') 
n=0 v 7 

% - Z v 


cos ncp 


cos ncp 


( 4 ) 


n=0 


rep 




n sin ncp 


n=0 


The terms in Q may be obtained from equations ( 4 ) by exchanging cosine 
for sine and vice versa and by making the right-hand side minus in the 
last equation for t . 


Displacements 


The radial and transverse strains 
from the stresses by Hooke f s law: 


e and e in the sheet follow 
r cp 
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\ 


E s € r - - a r - Vcfq) 


(S + s) 




( 5 ) 




J 


where E s is the modulus of elasticity of the sheet material. 

The stress-strain relations of the reinforcement ring can he expressed 
in the f orm ( see appendix A) : 


where, for the ring, e r and * r are, respectively, the extension strain 

and curvature, T is the normal stress resultant, G the tending moment 
(at any cp) , d the width of the ring, A^, and I r the cross-sectional 

area and moment of inertia of the ring (which may be functions of qp) , 
and u_k and u^, the circumferential and radial displacements of the ring 

at its outer circumference (r = aj_) . 

Since the displacements of the ring must be the same as those of the 
sheet where the ring is joined to the sheet, it follows that the displace- 
ments in equations (6a) and (6b) can be obtained by putting r = a^ in 
equations (5) for the displacements in the sheet. 

The stress and moment resultants T and G in the ring can be 
derived from the equations of equilibrium of the ring with the use of the 
fact that the unit loads acting on the reinforcement ring are due to the 
radial and shear stresses a r and t in the sheet along the circumference 

of the cutout ring. 


Cross-Sectional Properties in Terms of the Airy Function 

Substituting the expressions, as obtained in the preceding discussion 
(see appendix B) , for the displacements and the stress resultants into 
equations (6a) and (6b) and solving these equations for A r and I r yield 



(6a) 



( 6b) 


the following expressions for the cross-sectional area and moment of inertia of the cutout 
reinforcement ring: 


oo 



where c^, an( l c i ) 31,0 integration constants. Sine terms are included in equations (7a) 

and (7"b) by merely adding terms exactly similar to the corresponding cosine terms, with P n 
replaced by Q n . 


General Possibility of Exact Equivalence or Zero Constraints 

Equations (7a) and (7b) have the following significance. Suppose any type of original 
stress distribution as represented by an Airy function in the general form, equation (3a), is 
given (i.e*, given by the loading at the external edges of the sheet). This means that, for 
an equivalent cutout reinforcement, P n and are given. Then substitution for P n 

and Q n into equations (7a) and (7b) gives the distribution of cross-sectional cutout -ring 

area and moment of inertia required to give zero constraints, that is, required to replace 
with elastic equivalence the portion of the sheet cut out. If the values of A r and I r thus 
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obtained, are positive for any qp, then these are the correct and only 
values for zero constraints. If, on the other hand, the resulting values 
of A^. and I r are negative or infinite for certain values of cp, then 

an elastically equivalent reinforcement cannot be obtained. This procedure 
is considered in more detail in the following discussion. 


In equations (7a) and (7b), the P's (or Q's) in the general case 
represent the total stresses, that is, original plus constraint stresses, 
in the sheet. In case the reinforcement ring around the cutout is equivalent 
to the cutout portion of the sheet, then the constraint stresses will be 
zero and only the original stress distribution will exist in the sheet- 
Hence in that case, and in that case only, the P's will represent the 
original stresses. The original stresses in a sheet due to any loading of 
the external edges must be finite at all finite values of r including the 
center (r = 0) ; hence they can be represented in general by an Airy func- 
tion of the form given by equation (3a) in which the coefficients P n 
(or have the form (cf. equations (3 b)): 

\ 

P 0 = c 20 r 2 


P 


1 



> 


(3c) 


S n* 


= c 


In 


r 11 + c 


„n +2 


Rn 




The expressions (3c) can, if desired, be put into equations (7a) and (7b) to 
give slightly more explicit expressions for and I r (as functions 

of cp) in terms of the constants c which depend on the given edge-loading 
conditions. 


SPECIAL CASES 

Cases of Zero Constraint Stresses 


The problem of the choice of an exactly equivalent reinforcement ring 
around a circular cutout in a plane sheet is, in principle, completely 
solved in GENERAL EQUATIONS AND GENERAL SOLUTIONS. It remains to show some 
of the practical implications of this solution. For this purpose, the 
special case will first be considered in which the edge loading on the 
sheet is such that the original stress distribution can be represented by 
an Airy function with only a single trigonometric term, that is. 
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or 


ilr = P n cos ncp 


i = Q n sin ncp 


where P n and are given by expressions (3c) • This class of cases 

will lead to cross sections constant along the circumference of the ring. 


For n = 0, equations (7a) and (7b) 
radial symmetr^ lead to the result: 


(V 


= c-> = 0 'because of the 


Ar = 


at 


E c 


V E 


(8a) 


I 


r 



d _ a l 
5a atC2Q 



(8b) 


The fact that the constant c^ is arbitrary shows that in this case the 

moment of inertia I is arbitrary; but the cross-sectional area must, 
for zero constraints, have the value given by equation (8a). 

For n = 1 there results (assuming symmetry with respect to the 
y-axis and setting, therefore, C2 = 0) : 




where, for abbreviation, 
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In order that I r he finite, it is necessary that 


Whence , 



(9) 


Again, as for the case n = 0, I r may he chosen as desired, hut the 
area Ay required for zero constraints is fixed hy equation (9)» 

For n > 2, hy setting c q = c 2 = c 3 = 0 to obtain a constant cross 
section, equations (7a) and (7b) lead to: 


E s 

A r = at — 


Kn + n + 2 


K 


2a 


I r = a3 t -2 


E £ 

7 t 


o\ 



; n y (1+ v) 

+ (n + 1) [jn - 

2)v + n + 2 J + ^ jjn 

K 

1 

[j n " 1 ^2a " l 

] ♦ [ (n + 1] t - i 

n 2 - 1 

-n(l + v)k + 

[2 - n - V(2 + n)] 


(10a) 


(10b) 


where 


K = 


c ln 


c a 

3 n 1 


and K is prescribed hy the original stress distribution* Here both the 
moment of inertia and the area are fixed if zero constraints are to be 
achieved. 
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Two technically important cases of exactly equivalent reinforcements 
having constant cross sections are included in equations (8a), (8b), (10a), 
and (10b) and are those of hydrostatic stress and of pure shear. This can 
be readily seen in the following discussion. 

Homogeneous hydrostatic stress o r = a T r( p = 0- - In this case the 

Airy function, as may be verified by* * equations (2), is given by equations (3b) 
as 


* = P 0 = c 20 r " 


°r a cp ^ c 20 


From equations (7a) and (7b) it follows, as in equations (8a) and (8t>) , that 

(Ha) 


at E s 


Ay = 

r 1 - V E. 


I = Arbitrary value 


(Hb) 


These equations show that the equivalence of the reinforcement is assured for 
the definite value of the cross section Ap given by equation (lla) but for 

an arbitrary value of the moment of inertia I r of the ring. The value 

of I r will then, of course, be chosen as small as is compatible with* A p 
and with buckling considerations. 


Homogeneous shear distribution . - This (original) stress distribution 


can be expressed by 


T xy “ ^0 


in Cartesian coordinates. In polar 


coordinates. 


a p = -t q sin 2cp 

cr = t sin 2cp 
cp o ^ 

T = ”T_ COS 2 Cp 

rep o ^ 


( 12 ) 
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The corresponding Airy function, from equations (2) , is seen to he 

r 2 

^ = t q — sin 2cp (12a) 

signifying that (cf. equations (3c)) 



+ c 



where 


'12 


- T o / 2 


C32 _ 0 


Using equations (7a) and (7b), but with sine terms instead of cosine terms 
and with cj_ = C2 = = 0 to obtain a constant cross section, gives 


A r - 


at 




■‘-'(‘-a 1 ' 


(13) 


*r - 


- s) E fl 

6(1 + v) ®r 


(14) 


States of Stress with Imperfectly Equivalent Reinforcement Rings 

In the previous special cases of one trigonometric term for the Airy 
function corresponding to the original stress distribution in the sheet, it 
was seen that a reinforcement around a cutout could, theoretically, he 
designed so as to produce no additional stresses due to the cutout- The 
required cross sections were, in fact, constant along the circumference* 

For most other cases of external edge loading, however, it will he found 
that equations (ja) and (7b), intended to give full equivalence, will lead 
to physically impossible (i.e., negative or infinite) values of A^, and I r 

for at least some values of cp and that zero constraints are therefore 
impossible. Moreover, even In cases where zero constraints can be theore- 
tically achieved, it may be found (as in the case of pure shear, discussed 
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further in the following section) that the required geometric properties 
of the cross sections are in practice not realizable. In such cases it is 
of interest to determine how small the constraint stresses can be made for 
(preferably) constant cross sections of a reinforcement ring. 

The constraint stresses can be calculated without difficulty by 
substituting some given constant values of and I r into equations (7a) 

and (7b), using expressions (3b) for the P n (or Q n ) , including both the 
original (given) and constraint (unknown) terms, and then determining the 
unknown constants c^ n and c^ n (corresponding to the constraint stresses) 
so that equations (7a) and (7b) are identically satisfied for any angle qp. 

The constraint stresses will be calculated here under the assumption that 
they vanish at infinity. A number of constants then disappear in the 
expressions (3b)* This condition is selected to give practically negligible 
constraint stresses near the edges of the sheet. Because of this condition 
the constraint stresses in the sheet as derived in this report must be 
considered as approximate (they would be exact for an infinite sheet), but 
the approximation will be good if, as is commonly the case, the constraint 
stresses diminish sufficiently rapidly with distance from the cutout so as 
to have a very small or negligible magnitude at the edges of the sheet. 

Three special cases of edge loading are now considered in detail. 

Pure shear . - It will be found, upon closer examination of expressions ( 13 ) 
and (14), that the moment of inertia of the ring required for an 
exactly equivalent reinforcement in the case of pure shear is very high 
if its cross-sectional area be not higher than that required by equation ( 13 )* 
In actual design, therefore, the cross-sectional properties given by 
equations ( 13 ) and (l4) cannot ordinarily be realized, and constraint 
stresses must be allowed. In the following paragraphs the method of cal- 
culating these stresses is given in detail. 

For pure shear it suffices, except for the addition of terms in r, 
to employ only the trigonometric term in the series for the Airy stress 
function which was used in the case of no constraints, namely: 


\Jr = Q 2 sin 2qp 


where (cf. equations (3b)) 

^2 = ° l2 r2 + ° 22 r " 2 + °^2 

and, as before, 

_ T _o 
C 12 " 2 


(15) 


NACA TN No. 1852 


15 


The constraints are given by the terms in c 2 2 snd c k2t which are the 
unknowns. In expressions (3b), the constant c^ 2 must be made zero because 

otherwise it would give constraint stresses which do not decrease with 
distance from the cutout. Putting equation ( 15) into equations (7 a ) and (7b) 
and introducing for abbreviation 


-I 


. _ -4 , 

k 22 = c 2 2 a i > 


k 42 = V“l 


(15a) 


the expressions for A^. and I r become: 


A. = ^ 


Tq “ ^-^22 


r E r 1 + V 


t o (1 + 45) + (6 - 86)k 22 - (^ & 


(16a) 


V 1 + V, 


k 42 


^s a-’t 
r " E r 6(1 + v) 


3. T 0 ( 1 " 8) + (2 + 65)k 22 + (2 + 26)k 


'42 


T o " 2k 22 " ! + v k 42 


(16b) 


Let 


E. 


r ^r H 


at 


Eg T 2(1 + V) 


a 




a 3 t 


Eg r 24(1 + V)' 


Then solving equations (l6a) and (l6h) for k22 ^42 with v - 0*3 gives 


^22 


(1 - 5 - 0.251) (0-462 - 6. 156) a - (2 + 25 + 0-77i) [a(0. 5 + 2& ) " l] 

- T ... — ■ 

° (2 + 65 + 0.51) (6.156 - 0.462)a + (2 + 25 + 0.771) & + a(3 - 46)] 


(2 + 65 + 0.5i)[a(0.5 + 25) - l] - (l - 6 - 0.25l)[2 + a(3 - 46)] 
42 ° (2 + 65 + 0-51) (6.155 - 0.462)a + (2 + 26 + 0,77i)[.2 + a(3 - 45)] 




/ 


H 

o , \ 


r (17) 


Substituting the terms in c 2 2 and c h2 of equation (15) into the stress expressions and using the 
abbreviations of expressions (15a), the expressions for the constraint stresses become: 


^22 (a-,) + 4k 42 (a 


% = 6k 22(r^ 4 sin 2( p 


A 


sin 2cp 


> 


■<Pl 


*"(k 


-4 


2k 42( 

m 

cos 

2 9 

k 22 

and 

k 42 

can 


/ 


of equations (17), and the total stresses in the sheet can then be calculated by determining the 
constraint stresses as given by equations (l8) and adding these to the original stresses, 
equations (12). 


(18) 
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The following table gives the values of kg2 and- k^ 2 for several 

reasonably practical sets of values of thickness, area, and the stiffness 
ratios 6, a, and i. The entry values for 5=0 must be understood to 
be an approximation for a very small width of the ring. The approximation 
i = 0 for very small values of 5 introduces only negligible errors. 


PURE SHEAR 


a 

0 

1 

2 

4 

8 

5 = 0 ; i = 0 

’W’o 
*42 /t 

0.5 

-1.0 

0.161 

-.661 

0.0652 

-•565 



5 = 0 . 1 ; i = 0.1 

k W T o 

V/ T o 

— 

0.0505 

-.¥(-0 

-0.0682 

-.304 



5 = o.l; i = 0.2 

k 22/ T 0 

\s/ T o 

— 


-0.0675 

-.284 

- 0.154 

-.183 


6 = o.l; i = 0-4 

W T 0 
\ 2 l \ 

— 



- 0.153 

-.149 

- 0.226 

-.0884 
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It will be seen from this table and equations (l8) and (12) that in 
actual design it is difficult to avoid entirely some stress concentration 


at the cutout 



Without any reinforcement (i.e., for 5 = i = a = 


the stress concentration at the cutout is extremely high, for 


a 


cprp 


4-a 


*o 


where a is the total tangential normal stress with cutout and a is 

cp T ^ <P Q 

the original stress (without cutout) . It will be found that reinforcements 
with practical cross sections, such as those in the foregoing table, will 
relieve this normal-stress concentration, although they will introduce a 
smaller shear-stress concentration T rC p* From the point of view of minimum 

ratios of total stresses to local original stresses, the most suitable 
cross section given in the foregoing table is 5 = 0.1, i = 0.1, a = 1. 

The stresses for this section are given in the following table. (Also, 
see table I. ) 


PURE SHEAR 


5 = 0 . 1 ; Aj. = 


at _ a3t 

2(1 + V)’ r 240(1 + V)’ 


V = 0-3 


r/ a l 

%/ T 0 

v/ T ° 


Q 

Q 

0 ^ 

\l%o 

Tr<p r/ Trcp o 

1.0 

1.46 

0-303 

0-577 

-0.46 

1.303 

1-577 

1.25 

1.00 

.124 

.438 

.00 

1.124 

1.438 


Uniform single tension or compression . ~ In Cartesian coordinates, this 
state of original stress is given hy 



= a 


o 


o 


y 


o 



= 0 


where a is a constant, 
o 
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In polar coordinates, the same state of stress is: 


= -jp(2 - 2 cos 2cp) 



(19) 


T 



The Airy function following from equations (2) and (19) has the form 


It is easy to show that in this case, a perfectly equivalent reinforce- 
ment ring is not realizable (see appendix C). The Airy function for the 
final state of stress in the sheet will therefore contain terms representing 
stresses due to the reinforced cutout, that is, constraint terms, in addi- 
tion to those corresponding to the original state of stress. Referring to 
equation (20), it is seen that the Airy function for the original (super- 
script o) state of stress is of the form 





( 20 ) 


*0 = P 0° + P 2° 008 2< P 


where 



,2 



2 


c 


_ a o 

20 " c 12 = IT 


with 
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The constraint stresses, which will occur in the sheet with the reinforced 
cutout, are included by merely adding to the expressions for Pq° and P2° 

the additional terms appearing in expressions (3b). Thus, 


where 


= P 0 + P 2 008 


P 0 = c 20 r2 + c 30 lo Se r 
P 2 = c 12 r + c 22 r 2 + c 42 




( 21 ) 


c = c = - 
c 20 12 k 

in expressions (3b) can and must be set 

equal to zero, since otherwise they would give constraint stresses which 
do not diminish with increasing distance r from the cutout. The 
constant c^ has been omitted since it would not contribute anything to the 

stresses or displacements. In expressions (21), C2Q and c^2 are given 

/ oA 

(equal to — j by the original state of stress, while c^q, c 22> c 42 

are the unknowns. These unknowns are determined for a given reinforcement 
ring by putting expressions (21) into equations (7a) and (7b) , which express 
A^ and I r , and by satisfying equations (7a) and (7h) in every trigono- 
metric term. With the abbreviations: 


in which 


The constants Ci ir> and c „ 
40 32 
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c 30 a i " k 30 


c 22 a l 


-4 


- k 22 


c 42 a i" 2 = k 42 




|6(1 ♦ v) a ^ 

a-t 3 


Ay. E-v. 

at + 


\ 


J 


(22) 


the following result is obtained, (see appendix D) 


"22 


°o 

1 + aj_ 

^a-v) x 

2(1 + v) "2 

1.(4 ■ 1 + *k( 1 

0 

c 

+ 35 ♦ 1^(4 


1 [. R 2i 

— 1 + 5 + 


1 + V/ 


1 0^(1.+ 45) 



- 1 + 5 + 


21 

1 + V, 


j^a 1 (6 - 85) + 2j 


(23) 


k 42 = a o 


|(l + 38 + ij |l - <^(1 + 45)] - ^-(i-L -1+6) |a x (6 - 86) + 2 J 

C 1 + 38 * 1 ^( i£ rrf) - f 1 + 5 + rri)\y i6 - 88) + 2 ] 
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Putting expressions (21) into equations (4) and omitting the terms in c 2 
and c^ 2 , the expressions for the constraint stresses are seen to be: 



(24) 


In order to calculate the stresses in a plane sheet with a reinforced 
cutout under an external edge loading causing an original state of stress 
(i-e., state of stress if there were no cutout at all in the sheet) 
given by o x = o Q = Constant, it is necessary to calculate 6, i-, , 

and in accordance with expressions (22) from the dimensions of the 

reinforcement ring, determine k^ Q , k 22 , and k^ 2 directly from equa- 
tions (23), find the constraint -stress distribution from equations (24), 
and add this to the original stress distribution, equations (19). 


The case of uniform tension corresponds approximately to the experi- 
ments carried out for this report and, with the analytical calculation 
of the stresses in the test specimen (appendix G) , serves as an illustrative 
example • 


Pure -bending stress - - If an I-beam consisting of a rectangular web 
and flanges, as shown in figure 1, is subjected to a pure-bending moment M 
about an axis z perpendicular to the plane of the plate, then the stress 
distribution in the plate will be: 


a 



a 


y 


o 


x ^o 


= 0 


where 


NACA TN No. 1852 


23 


In polar coordinates this state of stress Is given by: 


- — r(cos 39 “ cos cp) 

O 4i w 


M \ 

r (3 cos cp + cos 39) 



> 


(25) 



It can he verified by equations (2) that the Airy function for this 
original state of stress is given by 


By substitution into equations (7a) and (7b) it is found (see appendix E) 
that it is not possible to design a reinforcement ring producing exactly 
zero constraint stresses. 

To determine the approximate additional stresses (in this case 
unavoidable) due to a cutout reinforced by a ring of given dimensions, 
there must, as in the preceding case, be added to the Airy function the 
additional terms in expressions (3b) corresponding to the same trigo- 
nometric orders as occur in Thus, for the stresses in the sheet 

under pure bending with a reinforced cutout 



where 


c 2i ' 3c 13 " “si. 


M 


W 


■f = Pj_ cos ip + cos 3qp 


P 1 = c 21 r3 + c 31 r 


( 26 ) 


P3 = + C2^ r ~3 + c 43 r 
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The constants and in expressions (3"b) have heen omitted since 

they would otherwise violate the assumed condition of constraint stresses 

C, 


decreasing to small values at the edges of the sheet. The constant 
has heen omitted since it has no effect on the calculations. 


"11 


Substituting equations (26) into equations (7a) and ( 7~b ) , determining 
the constants c^, c^, and c^ (c 2l and c^, as stated, are fixed 

by the applied bending moment) so that equations (7a) and (7b) are 
identically satisfied (see appendix F) , and using the abbreviations 


= 

a^t 


c 2 3 a l" 5 S k ' 


23 


c 4 3 a i" 3 = k ^3 


c 3i a i’ 3 S k 3l 


Maj_ 


2kl- 


= K 


W 


A r 

F b al 2 * 1 + v) 


'N 




= a 


? -r 2 Mi + v) s i , 
a 3t ' 


(27) 


and the assumption V =0.3* the following results are obtained: 

•(1 + &)(<& + 2k 31 ) 


k 3 = 


X 


k 4 3 = " K 


(i - 1 + 26) 

1 

[a(2 - 4.58) + lj 

+ (i + 1 + 46)[a(l + 4.58) - l 1 

I 

(l.36i + 1 + 26) 

[a(2 - 4.56) + lj 

1-1 

[a( 0.641 - 6-125) + O.333] 


1 + 46) 


k 23 = -K' 


(i - 1 + 25)[a(0.64l - 6.125) + 0- 33£] 

M 

[a(l + 4-58) - l] (l.36i + 1 + 28) 

(i + 1 + 45)[a(0.64l - 6.218) + O.333] 

l-l 

[a(2 - 4-58) + l] 

| (l.36i +1+25) 


(28) 


k u = -K 


Ja(2-78 + 40-65) - 3 + 65] + k2 3 [a(l8 - 40.45) + 9] + l% 3 [a(5.76 - 556) + 3] 

a (l " |) + 2 + 25 




Putting equations (26) into equations (4) and omitting the terms in ^21 &nd C 13 (which give the 
original stresses) , the constraint stresses in the sheet due to the reinforced cutout are found to he: 


q ’ ' a 31 ( q )’ 3 008 * ' [^ S ©' 


r\-3 


+ 10k 4 3 (^ 


COS 39 


= 2k Q1 (r) 3 COS CP + 


>j_ 3i 


|^ 2k 23 (^) 5 + 2k 4 3 f^)’ 3 J cos 3 * 


(29) 


rc P 1 ” 31 V a i. 


r\-3 


T — = 2k 01 fc) J sin cp + jl2k 23 Q 5 + 6k^ 3 1 ^— ) | sin 3cp 


S' 5 ] 


ro 

ui 
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With the constants k^, kg^, ^31 determined from equations (27) 

and (28) , the stresses in the sheet can be readily calculated by adding 
the constraint stresses, as given by equations (29), to the original 
stresses, as given by equations (25). 

In many actual cases, the reinforcement rings are quite narrow, so 
uhat if the approximations 5=0 and i = 0 are made only negligible 
errors will be introduced. The stresses in the sheet are then functions 
of only the cross-sectional area of the ring. This is true, of course, 
for any type of external edge loading on the sheet as well as for a pure- 
bending load. The following table gives the values of the constants k 
for pure bending for different values of the area concentrated in such a 
line reinforcement - 


PUKE BENDING 
[5 = 0; i = 0] 


K 

a = — 2(l + v) 
at 

0 

1/2 

1 

2 

V K 

3.00 

1-857 

1.48 

1.182 

k 2 3 / K 

-2.00 

1 

Oo 

ui 

-0 

-.486 

-.182 

2 W K 

-6.00 

-2.31 

.426 

3.22 


Figure 2 shows the stress distribution in the sheet for a = 1 and 
for comparison, for a = 0 (i.e., no reinforcement around the cutout). 


Summary of Special Cases 

Table I gives a brief summary of the numerical results of the 
theoretical investigation of the special cases treated here. In this 
table, W.E. denotes the ratio of the weight of material added (in order 
to form the reinforcement ring) to the weight of material cut out. Assuming 
that material of the same density is used for both ring and sheet, this 
weight ratio is given by 
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Table I shows also the ratios, at the cutout 




of maximum 


total (i.e«, original plus constraint) stresses to maximum original stresses- 
In the column headed ll Remarks, " the term "satisfactory" means that the 
cross section concerned has practical dimensions and that at the same 
time it prevents high stresses in the sheet. 


An examination of table I shows that it is theoretically possible to 
obtain, for the cases of hydrostatic stress, pure shear, single tension 
(or compression) , and pure bending, zero or minimum constraints and mi nlmum 
total stresses by reinforcement rings. However, except for hydrostatic 
stress, the required cross sections of such rings would have very high 
moments of inertia but very small widths and areas. Such sections are 
extremely difficult to design in practice. The table also shows, on the 
other hand, that in all cases considered here, not the minimum but at 
least fairly low total stresses in the sheet can be obtained by using 
appropriate reinforcement rings of practical dimensions. It appears, 
moreover, that in several instances, such rings may weigh less than the 
material removed from the sheet to form the cutout- 


TEST ON PLANE SHEET WITH REINFORCED CUTOUT 


For the experimental part of this research a test was made on a 
plane sheet with a reinforced circular cutout under a tensile load, 
realized by a heavy I-beam transmitting four concentrated loads to the 
sheet (see figs. 3 and. 4). The loads were produced by adjustable Jack- 
screws and measured by calibrated strain gages on the eight connecting 
links. 

Before carrying out the main test, a preliminary test was made with 
a sheet of the same dimensions but without cutout. The purpose of this 
preliminary test was to ascertain whether the means of application of 
the load would produce a uniform tension stress (a x = c Q J. Two different 

total-load stages were used, with the loads distributed as evenly as 
possible over the length of the I-beam. The total loads for the preliminary 
test were 1475 and 3016 pounds. For the main test (i-e-, with cutout), 
the loads were 3000 and 4020 pounds. 

The locations of the strain gages (including rosettes) together with 
the results of experimental measurements in the sheet without cutout are 
shown in figures 5 and 6* 

From the measured values of the test on the sheet without cutout, 
it is seen that the axial strain e x was distributed fairly uniformly 
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throughout the sheet. At the higher load, the maximum deviations from the 
average measured value were -10 and 6-4 percent while, from the expected 

value of € = = 2.98 X 10~\ the maximum 

^ tE s 0.052 x 18.5 X 10.5 X 10° 

deviation was 10.1+ percent. 

Near .the loaded edges rosette strain gages were attached to show the 
influence of the lateral constraint of the riveted Joint of sheet and 
I-beam. The smal 1 average magnitude of the transverse strain e v in 

the sheet near the loaded edges indicated that in fact the load was 
transmitted by the heavy I-hars in such a way that the transverse strain 

^and not the stress was practically zero at these edges. If the 

stress Oy had teen zero at the edges, then the magnitude of e would 

have teen much higher, namely e = -Ve = -0.3^ • 

y x x 

Having checked this type of original stress (or strain) distribution 
in the sheet obtained with the particular loading used here (see fig. 3), 
a circular hole was made in the center of the same sheet specimen and was 
reinforced by a ring of dimensions given in figure 3* Strain rosettes and 
simple strain gages were then placed on the sheet in the symmetric posi- 
tions indicated in figure 7* Two different loads were applied, 3000 and 
4020 pounds. Since both loads produced proportional results, only those 
corresponding to the higher load are given here. The results of the 
measurements are shown in figure 7* 

It may be of interest to compare the strains measured in the actual 
test specimen with the strains calculated, in accordance with the theory, 
on the basis of an infinite sheet. Since the boundary conditions of the 
test piece are obviously different from those of an infinite sheet, the 
strains may be expected to be different in these two cases. A comparison 
of the strains for these cases may nevertheless be instructive in indicating 
the influence which the finite-edge conditions have on the strain (or 
stress) distribution in the sheet. This comparison is given by the strain 
diagram in figure 7* Because of the symmetry about both x- and y-axes, 
the theoretical values for the infinite sheet are given only at the points 
indicated in the diagram. 

From figure 7* it ‘will be observed that, qualitatively, the finite 
sheet behaves quite similarly to the infinite sheet. For example, in 

both cases the transverse normal strain /€ \ decreases steadilv 

( 

a i ■ 
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•with 9 from a maximum at cp = 0° to a mi m'miim at cp = 90°. 

Also ( e x near the rigid "bars varies in both cases from a 

\ J Vx=9-2 in. 

maximum at the edge of the sheet to a minimum in the center. (See graph 
below strain diagram, fig* 7>) The radial strains were all of a relatively 
low order of magnitude, but it can be stated that here, also, both the 
order of magnitude and the variation of corresponding to the calcula- 

tions for an infinite sheet were similar to those of the test specimen. 


The effect of the particular boundary conditions of the test piece 
appears to be pronounced at two places. The maximum transverse normal 
strain (perpendicular to the radius and, at 9 = 0 , in the direction 

of the general tension stress) near the cutout ring = 1 . 38 ^ is higher 

for the test piece than for the infinite sheet, the percentage difference 


being 


5-0 6 - 4.07 

5.06 


x 100 - 19*6 percent. The axial strains near the loaded 


bars (see graph below strain diagram, fig. 7 ) decrease more sharply in the 
actual specimen than in the specimen calculated as an infinite sheet, the 


local percentage difference at the center being 


3.17 - 2.32 
2.32 


X 100 = 3 6* 6 per- 


cent, although the analytical values were practically equal to the experi- 
mental values at the quarter points of the sheet. The difference at the 
center must be due to the fact that in the actual test piece the rigid 
bars caused constant axial displacements and zero transverse displacements 
at the finite loaded edges, whereas in the theoretical work an infinite 
sheet was treated with constant axial stress at the remote loaded edges. 


The maximum strain measured in the sheet was the shearing strain 

at 9 = 45°, and it will be observed that this value ^ 7-66 x 10 was 
practically unaffected by the difference in the boundary conditions between 
the actual and the theoretical specimen, since the experimental and the 
theoretical values are practically equal. 


The calculations for the infinite-sheet specimen, based on the formulas 
developed in the analysis preceding the experiments, are shown in Appendix G. 

It may be remarked that in -these calculations account was taken of 
the rivet holes in the cutout ring. These holes had the effect of reducing 
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the actual cross-sectional area of the ring. To determine the amount of 
this reduction, a test with the same size and spacing of the rivet holes 
as in the sheet specimen was carried out. 

The ratio of the effective area to the full geometric area was thus 
found to he 0.80. Therefore in the calculation of the theoretical stresses 
in the test specimen used in the experiment the cross-sectional area of the 
cutout ring was taken as 0-80 times the full geometric area. 


CONCLUDING DISCUSSION 


The problem first considered was what must be the cross section of 
a ring reinforcing a circular cutout in a plane sheet in order that the 
stresses in the sheet remain unchanged by the cutout. The general solution 
to this problem is given by equations (7a) and (7b). In these equations 
the required distributions of moment of inertia I r and of cross-section 

area Ap along the circumference of the reinforcement ring are expressed 
in terms of the stress function (with coefficients P n and Q n as 

defined by equations ( 3 a) and ( 3 b)) for the original stress distribution, 
that is, for the stresses in the sheet without cutout or reinforcement 
rings. 

It was found that these expressions in terms of the circumferential 
angle cp for I r and lead to physically possible (i.e., positive) 

values only in a limited number of cases (for example, when the original 
stress function has no or only one trigonometric term, as in centric 
symmetry and in pure shear, respectively) . In the other cases, which 
include uniaxial tension and pure bending, constraint stresses, that is, 
additional stresses due to the reinforced cutout, are unavoidable. A 
method of calculating such stresses in the sheet for a given (constant) 
cross section of the reinforcement ring was developed, based on the 
requirement that the constraint stresses and displacements vanish suffi- 
ciently rapidly with increasing distance from the cutout so as to have a 
negligible influence at the edges of the sheet. The formulas, which are 
straightforward and convenient to apply, were derived in detail for the 
cases of pure shear, uniform axial tension (or compression), and pure 
bending. For example, in the case of uniform axial tension, it is 
necessary merely to calculate the values of the dimensionless constants iq 
and a-^ (proportional, respectively, to the cross-sectional moment of 

inertia and the area of the reinforcement ring) from the given data in 
accordance with the definitions, equations (22) , of iq and a^* The 

values of the constants k^o> ^ 22 > follow from the elastic 

properties of the ring and are determined by equations ( 23 ). The 
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constraint stresses (which must be added to the original stress distri- 
bution a y = o Q = Constant) then follow readily from equations (24) . 

In the formulas derived for the cases deemed technically important, 
it can he seen that for the narrow reinforcement rings commonly used only 
small errors will he introduced in the numerical calculations if the moment 
of inertia is taken as zero.l The results of the analysis therefore show 
that the stresses in a plane sheet with a cutout are a function only of 
the cross-sectional area of the ring reinforcing the cutout, while the 
moment of inertia is of practically no influence. This means that 
the ring experiences primarily tensile or compressive stress result- 
ants and not lending moments. 

In general the approximation given hy the formulas of the report will 
he closer the smaller the constraint stresses are at the edges of the finite 
sheet . 

In the experimental part, a test was made of a plane sheet with a 
reinforced cutout subjected to a uniform tensile displacement. The 
purpose of this experiment was to see how great an error is produced hy 
boundary conditions which differ from those for which the theoretical 
formulas are exact. In particular, the sheet specimen was (of course) 
not infinite, while at the loaded edges, the axial displacements, and not 
the stresses, were constant (with zero transverse displacements there) . 

On the two other opposite edges the normal and the tangential edge tractions 
were zero. 

It was found that qualitatively the strain (and therefore stress) 
distribution in the sheet specimen was quite similar to that predicted 
hy the theoretical formulas for an infinite specimen with the same rein- 
forced cutout. Quantitatively, the chief effect of the actual boundary 
conditions seemed to he at the loaded edges where the axial strains 
decreased more sharply toward the center of the sheet than predicted hy 
the formulas for the infinite sheet. 


The transverse strain 


°<p 


at the transverse center line of the 


sheet near the cutout 



was actually about 20 percent 


^This, of course, is true only if there is no concentration of the 
original stress along any radius of the cutout circle, a case which appears 
for a concentrated load at an external edge hut which was not considered 
in this report. 
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higher than according to the theory for the infinite sheet- In other 
respects, however, the quantitative results appeared not to he greatly 
affected by the particular boundary conditions of the test specimen, 
since the experimental values of the strains were fairly similar to the 
analytical values based on an infinite sheet. 


Polytechnic Institute 
Brooklyn, N. Y-, 


of Brooklyn 
June 27, 19^7 
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APPENDIX A 


STRESS-STRAIN RELATIONS OF REINFORCEMENT RING 


The theory of curved beams states that the change of curvature K r 
of a ring due to a bending moment G is 


K r = " 


E I 
r e 


where the effective moment of inertia I e is given by 


^■e r o 


br dr - Ar 


a i" d 


and where r = r Q denotes the radius of the neutral fiber, while b is 

the width perpendicular to the plane of bending of a cross section, and r 0 , 
moreover, is given by (cf., for example, reference l) 



The following table, however, shows that for a cross section of I-shape 
(for example) with values of d/a up to 0.4, a ring may be treated as a 
straight beam, so that it is permissible to put 

d 

r o = a l " 2 = a 


e 



I 

e 
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COMPARISON BETWEEN CURVED BEAM AND STRAIGHT BEAM 


d/a 

a l " | 
(in.) 

r o 

(in. ) 

*e 

(in. 4 ) 

X r 

(in. 1 * - ) 

0.2 

6 

5-973049 

0.041855 

0.04187 

•3 

6 

5-93044 

.119638 

.124966 

.4 

6 

5-865331 

.252761 

.260266 


It follows then that 


K r _ E I 
r r 


as in equation (6b). The extension strain € r and the curvature K r of 

a ring are given in terms of the radial and tangential displacements u~ 

r o 

and u_j_ at the neutral fiber r = r Q by the following well-known relations* 

s r a * % + “r D 


v 2 ■ -K * %). 


r=a n 


Since the radial displacement of the ring for any given cp will 

he constant along the width of the ring, it follows that 


"r" 2 * -C“r + %) r _ r - -(% * %•) 


The tangential strain, however, will vary along the width as, in fact, 
can he seen from figure 8, 
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Hence 


e a 
r 


■(** 


+ Uv. + U~ 


Thus, equations (6a) and (6b) follow. 


36 


NACA TN No. 1852 


APPENDIX B 

EXPLICIT EXPRESSIONS FOR STRESS RESULTANTS AND 
STRAINS OF REINFORCEMENT RING 


From figure 1 the following equilibrium conditions between the stress 
resultants T, G, and N in a ring and the unit load forces Y and Z 
can be derived: 


T + B - -Y^ - - (T. \ S] t (Bl) 

' / r=a 1 

B - T = -Z^ , (B2) 


G - Nr = -Ya,e = -fr \ a,te 
° 1 l rcp; r=ai 1 

d , , 

Noting that r © ~ a , e ^ (see appendix A), the last equation can be 
written as 


G - Na = - T^t 2 


(B3) 


Equations (Bl) to (B3) can he solved for T and G as follows: 
From equations (Bl) and (B2 ) , 


T + T = ai t(c r - T r(p ) r=£ 


(Bk) 


Moreover, from (Bl) and (B3)> 

G - -Ta - ai 2 t J(- r9 ) r=ai 


(B5) 


dqp + c± 
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By putting d r and t in terms of the stress function ^ (relations ( 2 )) 
equation (Bh) can he directly integrated with the results: 

for n 1, 


T = 


a^^t 


( r ' 1 *')r« 


+ c 


2 sin cp + c^ cos qp 


G = a l t | (^Mr^ 


a^ 2 t 


(r" 2 t) r=£ 


C 1 " c 2 a S ^ n 


c^a cos qp 


for n = 1, 


T 1 = 


C2 sin cp + c. 


cos qp 


Grj^ — “ 


sin cp - 


° 3 a 


cos cp + 


C 1 + 


,-l 


r 


-2 


1 * , )r=a 1 


Substitution for \(f hy means of relations (3) will give expressions 
for T and G in terms of P n and cp. 

Expressions in terms of P n for the extension strain and curvature 

can he obtained hy putting equations (2) into equations (5) and integrating 
the latter to find u r and u -^. 2 Thus, 


2 

Two arbitrary functions fj_(r) and f 2 ( qp) will appear as a result 

of this integration. However, it can he shown that, for compatibility 
between deformation values (involving also the shear strain 7 ) and 

the equilibrium conditions, the functions must have forms which are 
already included in the general expressions ( 3 &) and ( 3 h) for the Airy 
function. 
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E 


/ .. \ I f + 2? + ^ 

<> + v -J v 


dr + p (i|r + ^) - V(i' + \(r) ' 


E s( U r + *t + ”r 2i 


d\ d P + + t V /" d\ 

sj - saj — to + r* - ? i,v - 5jj 


V (| * 2 - a ♦)- - V 


By using the expression (3a) for i|r, the strains in terms of P n 
are obtained. 


and cp 
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APPENDIX C 


IMPOSSIBILITY OE zero constraints in case of uniform tension 


( a x = Constant^ 


From the expression (20) for \jr Q , 


P 


0 



P 2 



Substituting these values into (7a) and (7b), putting c 2 = C3 = 0 (for 

symmetry about both x- and y-axes), and simplifying give the following 
expressions for the distribution of cross-sectional area and moment of 
inertia: 

A = at — 1 + cos 2cp 

E r (1 - v) + (l + v + ^5(1 + v)] cos 2qp 


l r , a 3t h lil 5)(ltC ° S ^) 

®r 2(l - v) + 6(1 + v) cos 2qp 


where 6 = d/2a- 

It is not difficult to see that l r can be made positive for al i 

values of cp (by choosing a proper value for Cj_) but that A r neces- 
sarily becomes negative and infinite for some values of <p. Hence a 
perfectly equivalent reinforcement is not realizable in this case. 


4 =- 

O 


APPENDIX D 

DERIVATION OF EQUATIONS (23) PERTAINING TO AN ORIGINAL STATE OF UNIFORM TENSION 


Putting equations (21) into equations (7a) and (7b) with the abbreviations (22), the following 
equations result: 



i-i = 


K i + [if* 1 - 5 > + *22 * * 1 * + 6 (3 k 22 + ^s)] cos 2cp 


' 1 - V 
Ll2(l + V)°° 


k 30 

6 J 




22 


■ 1 ! v k «] 


cos 2cp 


where 

Kj_ = -k 1 a 1 1 + -^(1 - 5) + k 30 [(l + 5) log a^^ - &J 

Clearing each of the preceding two equations of fractions and equating the coefficients of equal 
trigonometric orders (viz, constant terms and terms with cos 2qp) , the following four relations 
are obtained: 
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2 + k 3 0 = a l 


p 0 (i -y) 

2(1 + V) *30 


2 — ^k 22 — ^1 


[°0 (I + 2& ) 


+ kgg( 6 - 8&) + k^ 2 


( 


K 1 = i 


1 - v k 30 

1 |_12(1 + v) a o " 6 J 


-(1-5) + k 22 (l + 35) + k 42 (l + 6) = i-,(r ■ k 22 




4v - l65 \ 

1 + v / 



The first of these four equations can he solved immediately for k^Q* The 
third equation determines the value of Kj_, which, however, has no influence 


on the stresses in the sheet. Hence this equation need not be considered 
any further here. The second and fourth equations can be solved simul- 
taneously for k 22 and k^g. The results are given by equations (23). 
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APPENDIX E 

PRACTICAL IMPOSSIBILITY OF PERFECT EQUIVALENCE IN 
CASE OF PURE BENDING 

From the expression for \|r Q , it is seen that, corresponding to the 
(original) state of stress due to pure tending. 


Substitution into equations (7a) and (7b) leads to the following expressions 
for Aj, and I r = Cg = 0, because of the antisymmetry about the x-axis 

and the symmetry about the y-axis): 






cos cp + cos 39 


at 


(2 + 95) (1 + v) 6 - 2V + 5(3V - l) 

(2 + 95) (1 + V) 


cos qp + cos 39 



/ 


(l - 26) ? cos cp cos 39 


r 24(1 + V) 


COS 39 


It will be shown now that and I r for zero constraints can both he 

made positive for any cp only if the value of the ratio 5 (identical 
with d/2a) is extremely low, too low in fact for practical purposes. 
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From the second of the two equations, in order that I r he everywhere 
positive it is necessary that 


a l 2t0 13 


= 6(1 + 5) 


Substituting, then, for into the expression for the area A r gives 


A? = 


at 


(2 + 98) (l + v) 


2(1 + 6) cos cp + cos 3qp 


6 - 2v + S(3V - 1) 
(1 + v)(2 + 98) 


cos qp + cos 3<P 


The value of this expression for Ap will he positive regardless of 
and only if 


2(1 + 5) 


6 - 2V + 5(3V - 1) 

(1 + V) (2 + 95) 


cp if 


With v = O.3, the positive root of this quadratic equation in 8 is 


6 = 0.00684 


In practical design, this low 
will he incompatible with the 

a3t 


I„ = 


24(1 + V) 


(1 - 28). 


value 

large 


for 8 (identical with d/2a, 
required moment of inertia 


fig. 
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DERIVATION OF EQUATIONS (28) PERTAINING TO CASE OF PURE BENDING 


Substitution of the expressions (26) for P-^ and P^ into equations (7a) and (7b), with the 
abbreviations (27) , leads to the following results • 

f-k 3 - (1 + 5) (6k + 2k 31 )J cos cp + [7(1 - 26) - k 23 (l + 45) - k^d + 28)J cos 3 <P 


1 = 


K + ^3 + suTTj^J 008 3<p 


k^ cos cp - |3K + 3k 23 + k^^J 

cos 3cp 

{■ 

’ 3(3 - v) | c „3 - 9v ' 
1 + v " 2(1 + V) 

+ k 3l( 1_ 1) 

J 

> cos qp + 

_-4 + ? 6 ) +1i 23 

( 27-N , k 2 + 1QV -98(5+ V) 

,(6--Sj +1 %3 2(1 + v) j 

cos 3T 


Clearing each of these two equations of fractions and equating coefficients of like trigonometric terms 
(cos cp and cos 3<p) give four simultaneous equations • 


-k 3 - (1 + 8)(6fe + 2k 31 ) = 


5 + v 


K(1 - 28) - kg 3 (l + 48) - k U3 (l + 28) = i | K + kg 3 + 3(l + ""v“ k 43 

k 3 -- a | K [- 'Trd + 8 2Trr?>J 

-( 3 K * 3k 23 * k 43 ) ■ - [-k( 3 ♦ - f ^ 


+ k 31 (l 


2 + 10V - 95(5 + v) 


2(1 + V) 
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The second and fourth equations are linear equations in 


‘23 


and 


V 


which can therefore he readily solved for these two unknowns. By 
eliminating k^ from the first and third equations, k^ can be solved 

for in terms of k 0 -, and 


23 


k 4V The results are given by equations (28). 
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APPENDIX G 

CALCULATION OF STRAINS IN INFINITE -SHEET SPECIMEN WITH SAME 
REINFORCED CUTOUT AS ACTUAL SPECIMEN 


The dimensions of the specimen are given in figure 3- From the 
dimensions, 

a = 3. 66 inches 

a, =3.5+ -“ = 3* inches 
1 16 

t = 0-052 inch 

5 = — = — = 0-0426 

2a 2 x 3.66 


0.80 



X 0.434 


= O.IO85 square inch 


M 1 + v ) = 0.1085 x 1.3 = 0>t4 
at 3.66 x 0.052 


= 


f O x 0.434 
12 


= 1.1 X 10 ^ inch 14 ' 


i x = 6(1 + V) 


I r _ 7.8 x 1.1 x 10~ 3 
a^t (3.66)^ x 0.052 


= 3.37 x 10 


From equations (23) substitution of 
following numerical results for the 


the preceding values leads to the 
constants k^ 0 , k 22 , and k^: 
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k J0 - -0.173c o 

lt £2 = 0.0228o o 

\ 2 - -°- 2 * * * 52»o 


Frcan equations (24) it then follows that the calculated constraint 
stresses in the sheet will he: 



The calculated stresses without cutout (see equations (19)) are: 


a = — (l - cos 2cp) 


• o 2 


2 ^ + cos 


T rcp = ~ 3in 2c P 


Hence, adding the constraint to the original stresses, the total stresses 

in the sheet at any point (r, <p) are found to he: 


CT r = a o < 


0.500 -0.173(q)' 2 - lo. 500+0.1368^^ 





) - !• 00 8 (£ 

0 J 

cos 2qp 


°cp ■ M °-50° + °- 1 ' ? 3(4)' 2 + 


0 • 500 4 - 0 • 




-4 


cos 2cp, 


rep 


= Or 


r _x-2 
1/ 


0-500 + 0.504(^j- J - o*1368(— 


r_x"4 
1 / 


sin 2cp 
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From these stresses, the strains expressed in polar coordinates at any 
point are: 



_ 2(1 + V) 

7 rcp ~ Eg T rcp 


The strains e x in the direction of the load were obtained by first 
determining the Cartesian stresses o and a according to the relations: 


1 

2 


(°r + ' ( a r " °cp) 008 2( P + ^rcp sln 2 ? 


1 

°y 2 


°r + a qv) + ( a r “ a cp) 008 2c P ' 2r rcp sin 2c P 


The strain e T is then readily obtained by 


- £(°X - %) 


The value of the original stress o Q without cutout is 

Iff 


°° 2b t 


where 2b is the width of the sheet. For the test specimen and for a 
load of 4020 pounds, 

4020 

°° ~ 18.5 x 0.052 


= 4160 pounds per square inch 
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For the strain this can he expressed as: 


4l60 

10-5 x 10 6 x 4-02 


= 0.99 x 10 


-4 


The actual measurement gave 1.04 x 10~\ that is, an average inaccuracy 
of 5 percent in the strain gages. 

The numerical results of the preceding calculations are given in the 
data diagram (fig- 7) • 
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TABLE I 

SUMMARY OF THEORETICAL INVESTIGATION OF SPECIAL CASES 


Loading 

_d 

2a 

at 

* 

a 3 t 

W.R. 

(1) 

Ratios of maximum 
total stresses to 
maximum original 
stresses at r/a^ = 1 

Remarks 

( 3 ) 

a r 

0 

JL 

0<p o 

T rqp 

Tr< p 0 

°*o 

(2) 

Hydrostatic 
a r = Constant 

0 

1.13 

Any 

2.86 

1.00 

1.00 

1.00 

2 

Zero constraints; 
satisfactory 

.1 

1.43 

Any 

3-04 

1.00 

1.00 

1.00 

.2 

1.43 

Any 

3-21 

1.00 

1.00 

1.00 

Pure shear 

t = Constant 

.2 

.428 

0.1025 

.0875 

1.0 

1.0 

1.0 

1 

Zero constraints, but I r 

too high to conform with 
Aj, and d/a 

.1 

.385 

.0032 

.457 

-.16 

1.30 

1.58 

Satisfactory 

.1 


•1155 

.865 

1.0 

1.0 

1.0 

Zero constraints, but I r 

too high to conform with 
A^ and d/a 

Single uniform 
tensile ( or 
compressive) 
stress 

cr x = Constant 

.1 

.882 

•115 

1.69 

.822 

.872 

1.31 

3 

Low constraints, but I 

too high to conform with 
A r and d/a 

.2 

•795 

.1025 

1.232 

.783 

.790 

1.43 

Low constraints, but I 
too high to conform with 
A r and d/a 

0 

.615 

0 

1.23 

• 750 

1.11 

1.63 

Satisfactory 

.1 

.161 

0 

.682 

.512 

1.31 

1-59 


Satisfactory 

Pure bending 

0 

.385 

.0320 

• 770 ^ 

1.015 

.995 

1.015 

2 

Almost zero constraints, but 
I r too high to conform 

with A r and d/a 


0 

.385 

0 

•770 

1.01 

1.10 

1-34 


Satisfactory 


W-R. ratio of the weight of material added (to form reinforcement ring) to the weight 
of material cut out. 

2 

No reinforcement. 

Satisfactory means that the cross section has practical dimensions and that at the 
same time it prevents high stresses in the sheet. "Zero constraints" means that the rein- 
forcement ring is exactly equivalent to the portion of the sheet which has been cut out. 


^\NACA^^ = ’ 



VJ1 

ro 


Figure 1.- Stress resultants and displacements in a reinforcement ring represented in a 

polar -coordinate system. 


I 

> 

O 

H 

00 
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Figure 2. - Total-stress distribution in sheet produced by pure bending. 


I r = 0. m = 


°r 1 + a r 


0 _ + *9 o T r<P! + T 

- n = — , p 


r< P. 


max 


<T> 


max 


rep 


£. - = A - n 

’ 5 2a °* 


max 
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Figure 3 .- Experimental setup. 



Figure 4.- Test specimen, 
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369 

746 


364 

764 


360 

729 


369 

746 


9.2 


9.2 


A i 

i i 

1 h 

1 *7 ** 

III M 

liM 

HB 


AAA A 

1.56 1.52 1.56 1.59 

3.12 2.98 2.99 3.12 


1 3" 

>4 l6 


n n n n ^ 


20 . 6 " 


1.56 1.73 1.77 1.60 

3.24 3.29 3. 18 


I .48 
2.91 


1.61 1.59 1.65 1.68 

3.20 3.02 3.06 3.17 

A AAA 


363 

379 

760 

770 

Measured strains 

€ x 


values are for a total load 4P of 1475 pounds; bottom values are for 
a load of 3016 pounds. Rosette, A ; simple strain gage, p r | . Drawing 
not exactly to scale. 
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P P 


P P 


A 

.A 

A 

A 

- 0.09 

- 0.07 

- 0.08 

- 0.26 

-.20 

-. 15 

-.15 

-.60 


- 0.30 

- 0.08 

- 0.05 

- 0.08 

-.60 

-. 14 

-.13 

-.17 

A 

A 

A 

» 

A 



P P 


P P 


4 

Figure 6.- Measured strains x 10 transverse to direction of loading. 

Top values are for a total load 4P of 1475 pounds; bottom values are 
for a load of 3016 pounds. Rosette, a . For dimensions see figure 5. 
Drawing not exactly to scale. 
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Figure 7.- Strain measurements on plane sheet with reinforced cutout. 
Total load, 4020 pounds; rosette, A ; simple strain gage, f^l J values 
in parentheses are corresponding theoretical values for reinforced 
cutout but for an infinite sheet (see appendix G). For each rosette, 
the top, middle, and bottom values are, respectively, normal transverse, 

normal radial, and shear strain times 10^. For each simple gage, 
values given are those of axial strains e x 10^. 

A 
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Figure 8 .- Contribution of curvature of cutout ring to extension at outer 

boundary of ring. u^. = ^ “ ^r~ • 

-L o3. 
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